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Abstract. 

In view of the current efforts to extend the KASCADE experiment (KASCADE- 
Grande) for observations of Extensive Air Showers (EAS) of primary energies up to 
1-10 18 eV, the features of muon arrival time distributions and their correlations with 
other observable EAS quantities have been scrutinised on basis of high-energy EAS, 
simulated with the Monte Carlo code CORSIKA and using in general the QGSJET 
model as generator. Methodically various correlations of adequately defined arrival 
time parameters with other EAS parameters have been investigated by invoking 
non-parametric methods for the analysis of multivariate distributions, studying the 
classification and misclassification probabilities of various observable sets. It turns out 
that adding the arrival time information and the multiplicity of muons spanning the 
observed time distributions has distinct effects improving the mass discrimination. A 
further outcome of the studies is the feature that for the considered ranges of primary 
energies and of distances from the shower axis the discrimination power of global arrival 
time distributions referring to the arrival time of the shower core is only marginally 
enhanced as compared to local distributions referring to the arrival of the locally first 
muon. 
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1. Introduction 



Since the first experimental studies in 1953 by Bassi, Clark and Rossi and Jelley and 
Whitehouse arrival time distributions of the charge particle components of Extensive 
Air Showers (EAS) have been often theoretically and experimentally investigated under 
various aspects [U, H U, 0, H § ■ Systematic studies are performed particularly at 
large EAS detector installations: in Haverah Park 0, |ll], 0, 0, HI , at the 
Potchefstroom University H, in Akeno and on Mt. Chacaltaya Jl8|,[l9|,|0], at MSU |2l[| , 
more recently with the GREX/ COVER-PLASTEX [||,||,||,|5|,|| and KASCADE 
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q1\, |32|, |33], [34| experiments. Measurements of arrival time distributions of 



high energy hadrons near the shower core have been reported by the Maryland University 
group |35[] and by KASCADE |]36| . In addition to the fact that the observed arrival 
time distributions of the EAS particles display the phenomenological appearance (shape 
and structure) of the EAS disc, they provide also a coded picture of the longitudinal 
EAS development. Especially the muon component, when multiple scattering in the 
atmosphere at higher muon energies gets negligible, maps the distributions of the 
production heights via the time-of-flight of the muons from the loci of decay of the 
charged parent pions to the detector [27], [57]. Fig. 1 sketches the relation in a simplified 
way. Muons released in higher atmospheric altitudes (and observed on ground at fixed 
larger distances from the shower axis) show smaller delays relative to the arrival of 
the shower centre ("light front"), and their (relative) arrival times let expect some 
discrimination features for the mass of the EAS primary [^,[|[7J,^7J,^J. Thus the arrival 
time distributions of muons, produced in iron induced showers e.g. should get shifted 
to shorter delays due to the faster development as compared to proton induced showers 
of the same primary energy. Of course, a serious analysis of this mapping has to invoke 
realistic simulations of the EAS development and of the muon tracking, taking into 



account the influence of multiple scattering and off-axis production |37| of the muons. 
Fig. 1 indicates also that observations of the angle distributions of the muon incidence 
provide alternative experimental possibilities |38|,|39|], in case that multiple scattering 



effects do not seriously obscure that information [|37j]. The information content of the 
combination of both types of EAS observations, known as "Time- Track- Complimentary- 



Principle" [j3~7, IQ|, has been scrutinised in Ref. [29]. 



Arrival times t}, < tJ, 



< Tf. 



of the EAS muons, locally registered by timing 



detectors at a distance from the shower axis have to refer to a well defined zero- 



time, usually the arrival time r c of the shower core (global arrival times): 



glob 



e.g. 



Often there are experimental difficulties to determine the arrival time r c with 
sufficient experimental precision, and only "local" times are considered. They refer 
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Figure 1. Sketch of a simplified relation between muon arrival time and angle of muon 
incidence, respectively, with the production height, neglecting multiple scattering and 
off-axis production (/3 M c = velocity of the muon). Additionally (right) the definitions 
of the arrival time quantities characterising the shape and the structure of the muon 
disc are indicated. 



to the foremost (first locally registered) muon: 

Local arrival time distributions display the internal structure of the shower disc 
(characterised by various time parameters), but they do not carry information about 
the shape (curvature) of the front. Nevertheless studies based on simulations of the EAS 
development |16|,^9| have shown that mass discrimination effects are just pronounced 
in observations of At}, glob . 

For event-by-event observations with a fluctuating number of muons, the single 
relative arrival time distributions can be characterised by the mean values Ar mean , and 
by various quantiles Ar q , like the median Ar .5o, the first quartile Aro.25 and the third 
quartile Aro.75. For ordered statistics of measured times At\ < Ar 2 < ... < Ar„ 
and k := na + £, k integer and £ £ [0,1), the a-quantile Ar a (a £ (0,1)) is the 



Ar„ 



following pH2£ 

' (Ar fc + Ar k+1 )/2 : for £ = 
v Ar fc : for £ £ (0, 1) 

That means: in the case of large n, a fraction a of muons has arrival times less than 
Ar a . The mean values and dispersions (standard deviations) represent the time profile 
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of the EAS muon component. Measurements of muon arrival time distributions and the 
determination of the distributions of various different time quantities have been a subject 
of the current investigations of the KASCADE experiment [f2^,|2^,^0|,^T],|3^,^3|,^3|. In 
these investigations the temporal EAS structure has been studied in detail in dependence 
on the shower size N e , the muon number N^, or the truncated muon number Nff (which 
is used as approximate energy identifier in KASCADE observations jnj), respectively, 
and from the angle of EAS incidence 9. Special features arising from the observed 
muon multiplicity have been revealed [031 . More recently [34] using larger samples of 
shower events, experimentally observed with KASCADE, the sensitivity of local muon 
arrival time distributions and of their correlations with other EAS observable to the 
mass of the EAS primaries has been investigated. Methodically advanced statistical 
analysing techniques have been applied [fHql , based on Bayesian decision making for 
the classification of the EAS events. The procedure requires simulated distributions 
(processed through the detector response) as reference patterns, provided by the Monte 
Carlo simulation program CORSIKA j|3) with a particular model of the hadronic 
interactions. Thus the analysis necessarily implies some model dependence. But the 
consistency with the invoked QGSJET model J44| could be established thanks to the 
measurements with varying distance from the shower centre and varying muon 
multiplicity thresholds for being accepted in the observation sample. Nevertheless, it 
has been also shown that the mass discrimination power of local muon arrival time 
distributions is rather marginal in the observed range < 100 m and primary energies 
Eo < 10 16 eV, and it does not significantly help for the classification. 

The phenomenological features of the time structure of high energy EAS at larger 
distances from the shower axis have been experimentally studied with the Haverah 
Park |g,|g[g,|rj,|rj,|r5|,|T§ and the Akeno air shower arrays fg^g, and the 
results have been compared with phenomenological model predictions (scaling models 
and multiplicity prescriptions for the particle production in hadronic collisions). An 
analysis of the mass discrimination power on basis of modern models of the hadronic 
interactions has not yet been performed for these cases. 

The recently started extension of the KASCADE detector installation: KASCADE- 
Grande ||5| will not only provide the possibility to extend the studies of muon arrival 
time distributions and their correlations to larger distances R^ and larger energies of the 
observed EAS. There is also a chance to measure global arrival times, at least for R^ up 
to 350 m |46| . The investigations of this paper use the analysis techniques of our previous 
studies [26 , 34] arid refer mainly to the QGSJET [44|, in some few comparisons also to 
other models. The studies are focused to explore the basic information content of muon 
arrival time distributions of EAS with energies up to ca. 1-10 18 eV. The analyses are 
completely based on data of simulated showers and consider ideal cases, i.e. neglecting 
the influence of the detection system. The view of interest is the discrimination of the 
mass of the primary particles (with adopting a particular hadronic model as generator 
of the reference patterns). For some observable combinations the differences in the 
predictions of Monte Carlo simulations using different high-energy models as generators 
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have been scrutinised. Only marginal differences for the arrival time distributions have 
been found. 

2. Time profiles of the EAS muon component 

The present studies consider predictions of realistic and detailed Monte Carlo 
simulations of the EAS development initiated by cosmic ray primaries of different mass 
in the primary energy range up to 1-10 18 eV. The simulations have been performed by 
use of the program CORSIKA (version 5.64) ||43|| , invoking in general different models 
for the hadronic interaction: QGSJET (version 98) @, VENUS (version 4.125) @ 
and SIBYLL (version 1.6) jy|. For the particle interactions below E^ = 80 GeV 
the GHEISHA [^9| option is used. Earth magnetic field, observation level and particle 
detection thresholds have been chosen in accordance with the observation conditions of 
the KASCADE experiment, but without a detailed account for the detector response. 
The U.S. standard atmosphere (see Ref. has been adopted for the simulations. 



For a realistic description of the electron-photon component, instead of the NKG 



approximation |5(| the EGS option [51] has been preferred. 

A first set (I) of simulations comprises samples (of ca. 500 events each) of proton and 
Fe induced EAS of vertical incidence for various primary energy ranges, 10 15 eV, (1.78- 
3.16)-10 16 eV and 3.16-10 17 eV, calculated with the CORSIKA version 5.62. Another set 
(II) of samples, calculated with a spectral index of the power-law slope of —2.0 for eight 
energy ranges ( MO 16 - 1.78-10 16 eV; 1.78-10 16 - 3.16-10 16 eV; 3.16-10 16 - 5.62-10 16 eV; 
5.62-10 16 - MO 17 eV; MO 17 - 1.78-10 17 eV; 1.78-10 17 - 3.16-10 17 eV; 3.16-10 17 - 5.62-10 17 
eV; 5.62-10 17 - MO 18 eV) with the CORSIKA version 5.64 and the QGSJET model, 
comprises smaller number of events (100 or up to only 9 for the highest energy range), 
distributed randomly over an angle-of-incidence range of 0° — 30°. 

For each simulated EAS event a number of EAS observable have been reconstructed: 
the electromagnetic shower size N e , the muon content Nj[ , the shower age s, and in 
particular the time parameters Ar q of the muon arrival time distributions and their 
variations with the distance from the shower axis. For higher primary energies 
and approximating the experimental conditions of KASCADE-Grande, instead of N e 
the quantity N ch (the total number of charged particles) is preferentially considered. 
Consequently we deduce also an age quantity s c h derived from the lateral distribution 
of the charged particles. As approximate energy identifier instead of Nj[ the quantity 
Pti {Rfi — 600 m) with E t hr = 240 MeV is introduced, which has been often considered 



in the past for this purpose |52"|, |53[| 



The muon arrival time distributions have been calculated for muons with an energy 
threshold Ethr — 2.4 GeV, observed with a multiplicity n > 3 per event. The value of this 
energy threshold is chosen according to the muon detection threshold of the KASCADE 
Central Detector which is foreseen to be used for the time measurements ||34|| . Just for 



an impression average arrival time distributions for proton and iron induced EAS at 
two different energies and different muon energy detection thresholds are displayed in 
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Fig. 2. The figure indicates that differences for different kinds of primaries are, if ever, 
obvious in the initial part of the distributions, represented by the first quartile Ar .25- It 
gives also an impression about the order of the needed time resolution (about 1.5 ns) to 
reveal the differences. But generally the present paper is not focussed to the discussion 
of instrumental effects. 

The distributions of the different quantiles Ar q of event-by-event observations have 
been shown [^, |3]J to follow fairly well a phenomenological parameterisation by Im- 
probability distribution functions, with parameters varying with the primary energy 
and the distance R^. 

The following figures (Figs. 3 and 5) display the time profiles of the muon 
component, i.e. the variation of the mean < Ar q > of the time parameters and of their 
fluctuations (characterised by the standard deviation a q ) with R^. They are calculated 
with samples of set I and analysed for radial ranges up to 300 m (in bins Ai? M = 10 m). 




5 10 15 20 25 30 35 40 45 50 



Ax 9lob [ns] 




Figure 2. Average global muon arrival time distributions to be observed at the 
distances i? M to the shower core of proton and iron induced EAS at two different 
energies E = 10 15 eV (vertical incidence) and E = (3.16-5.62)-10 17 eV (0° -30°) and 
with the muon energy detection threshold E t h r —240 MeV (according to the threshold 
of the KASCADE detector array fl|). 
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Figure 3. (a) Comparison of local and global time profiles < Aro.50 > of the 
EAS muon component, induced by proton and Fe primaries of the energy range of 
(1.78-3.16)-10 16 eV, based on EAS Monte Carlo simulations using CORSIKA and the 
QGSJET model as generator. 



Fig. 3a compares local and global time profiles < Ar 50 > with respect of 
differences of the primary mass, while Fig. 3b displays the time profiles of different 
quartiles of the proton case. Obviously the local and global profiles do not exhibit 
pronounced differences for different primaries. The quantities < Ar 9lob > and < At 1oc > 
differ by an offset (Arf' ofe ) which varies with R^. It has been shown that differences of 
< Arf lob > are in the order of 1 ns for different primary muons in the considered range 
of distances from the shower axis |3| . 

Since the lateral distribution of the muon component displays also some differences 
for different kinds of primaries, the profile of a combined quantity ("reduced time") 
^ r o.50 = A r o.5o(-R/J/p^(-R^) is °f interest where p* represents the density of muons 
with energies > 2.4 GeV. This quantity is measurable with the KASCADE Central 
Detector and is related to the multiplicity of muons observed in the EAS event. 

As an example in Fig. 4 the fluctuations of p* are displayed for different distances 
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Figure 3. (6) Local and global time profiles of different quartiles < Ar g > of the EAS 
muon component, induced by protons of the energy Eq = (1.78-3.16)-10 16 eV from 
EAS Monte Carlo simulations using CORSIKA and the QGSJET model as generator. 



from the shower core, indicating an increasing separation of the distributions with 
increasing R^. This feature is reflected in the distributions of the time quantity Atq ^ ob . 

Fig. 5 compares the local < AtqJ™ > and global < Atq ^ oh > profiles. A 
marginally slight improvement of the discrimination of the global quantity for distances 
larger than 100 m is indicated. This feature gets more pronounced at higher energies. 
Though global quantities include the curvature of the shower disc, while local quantities 
display only the internal structure of the EAS. At not too large distances from the 
shower axis the offset between global and local profiles is nearly mass-independent. 

Generally we note that differences of distributions arising from different mass 
primaries are in the order of nanoseconds, increasing with the distance from the shower 
axis. For the energies < (1.78-3. 16)-10 16 eV (set I) also some exploratory calculations, 
adopting the VENUS or SYBILL model, have been performed and show that the model 
dependence leads to differences of the distributions (within the considered range) in 
the order of 2% with no clear trend [55 |. 
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Figure 4. Fluctuations of p* (left) and the distribution of the reduced time Atq 5 9 ° 
(right) of proton and iron induced EAS (Eq = (1.78-3.16)-10 16 eV) at different distances 
of the shower core. 



3. Correlated distributions 

The concept of modern cosmic ray experiments like KASCADE or KASCADE-Grande 
with a multi-component detector array is to deduce the information from correlated 
measurements of a larger number of observable EAS parameters for each individual 



event (see Ref. [pq] ). Particular observable quantities are the electromagnetic shower 



size N e , the muon content iV M (or the truncated muon number Nj[ in case of KASCADE) 
and other specific observable quantities, characterising the various EAS components. 
In view of the measuring possibilities of the KASCADE-Grande array the KASCADE 
observable N e , and iV* r (the latter loses also the role as approximate energy identifier 
for larger energies) are not accessible. For KASCADE-Grande they will be replaced by 
the total number of the charged particles N c h and a (partial) muon number, as measured 
with the original KASCADE detector array being embedded in the KASCADE-Grande 
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Figure 5. Reduced local and global time profiles < Atq J§ c > and < Ar 5 q > of 
the EAS muon component induced by proton and Fe primaries of the energy range 
of (1.78-3.16)-10 16 eV based on Monte Carlo simulations using the QGSJET model as 
generator. 



installation. This muon number is dependent on the distance to the shower core. In 
the present studies, for the sake of simplicity, we represent it by p^(i? M ) i.e. the muon 
density (with E thr = 240 MeV) at i? M and consider in particular p^(600 m). 

The variation of the muon density p M (i?^) with the primary energy is illustrated in 
Fig. 6 by results of the present simulations. It is obvious that p M (600 m) as well as the 
muon density at other distances (see Fig. 6) do not provide a strictly mass independent 
energy identifier. In addition as known for Nj[ , there are considerable fluctuations of 
p M (600 m) at fixed energies. 

The current and foreseen analyses of the experiments attempt to infer from the 
correlations of various observable an estimate simultaneously for the energy and the 
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Figure 6. The variation of the primary energy Eq with the average of p M (600 m) 
/9 M (300 m), ^(180 m) and ^(90 m) for proton, C and Fe induced EAS. 



mass of the primary cosmic particle inducing the EAS event. As the most powerful 
correlation the N e — iV* r correlation has been proven (see Ref. [|56|). It is expected that 
the N e — p At (i? /J ) correlation will exhibit a similar discrimination power which can be 
defined in terms of the true classification probabilities resulting from the analysis (see 
section 4). The features are indicated in Fig. 7 which compares the N e — Nj[ correlation 
with the N c h — p^(600 m) correlation at two different primary energies. 

Correlations with additional observable EAS parameters, though of relatively small 
own classification power, shrink the influence of the natural EAS fluctuations and test 
the consistency of the analysis |57|] . 

The observation of muon arrival time distributions at particular distances from 
the shower centre is combined with the observation of the local muon density p*(i? M ) 
with E thr = 2.4 GeV reflected by the muon multiplicity i.e. the registered number 
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at different primary energies. 



of muons spanning the single arrival time distribution. It turned out |)5] that the 
Ar g (i? M ) — p* {Rfj) correlation improves the mass discrimination and can be fairly well 
replaced by a combined parameter Ar*(i? M ) = Ar g (i2 M )/p*(i2 At ). This result has been 
already anticipated in the presentation of the EAS time profiles in section 2. 

Fig. 8 displays the correlation distributions of AtqJq c (R ij ) and Atq 5 Q Zob (i? M ), 
respectively, with the log 10 N e for proton and iron induced showers of the energy 
Eq = (1.78 — 3.16) ■ 10 16 eV in two different ranges. 

Special interest arises from the correlation with the so-called shower age s. This 
particular EAS parameter associated with the soft EAS component is supposed to carry 
also some information about the longitudinal EAS development. In the cascade theory 



(NKG approximation) |49] s is related to the actual atmospheric depth of the shower 



development via the lateral distribution involving the Moliere radius R m oiiere (which is a 
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Figure 8. Correlation of the local quantity Atq §a°{Rp) and of the global quantity 
Atq 5Q° b (R^) with the logm N e for proton and iron induced showers of the energy 
E Q = (1.00 - 1.78) • 10 16 eV in two different i? p ranges. 



characteristic unit of length of the scattering theory) and the ratio E /eo of the incident 
energy to the critical energy. However, the use of the age parameter resulting from 
the simplified handling of the electromagnetic component by the cascade theory does 
not describe realistically the electromagnetic component. Rather a full Monte Carlo 
simulation of the electromagnetic component (by use of the EGS option |jT] in the 
CORSIKA code) calculating the total intensity distribution of the electron component 
and the lateral electron distribution with all fluctuations, has to be required in order 
to obtain results being comparable with the reality. In order to derive from lateral 
distributions resulting from the EGS - Monte Carlo simulation an "age parameter" value 
(which is not defined within the Monte Carlo approach), the resulting lateral distribution 
is subsequently fitted by the NKG function /(i?/i? mo ;^ re ; N e , s) in the same manner as 
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Figure 9. Correlation of Atq ^° b {R^) with the age s and age s c h, respectively, for 
proton and iron induced showers of the energy Eq — (1.78 — 3.16) • 10 16 eV and 
E n = (1.78 — 3.16) • 10 17 , respectively, in two different i? M ranges. 



the experimental observations are analysed (R mo iiere = 78 m). Thus a (lateral) "age" 
is extracted ||58[| . In a similar way the lateral distributions of the charged particles are 
processed for deducing the age parameter s c h- However it should be noted that such 
a procedure is considered to be only a first approximation. The lateral distribution of 
charged EAS particles at larger distances from the shower core deviates from the NKG 
function, and an improved parametrisation should be used |59[| . Nevertheless Fig. 9 
displays the correlation with the age parameters, defined by the NKG function. They 
emphasise the suggestion to exploit in analyses of experimental data correlations with 
the age parameters, which are obviously related to the longitudinal development and 
could lead to an improved discrimination of protons and iron induced EAS. 

In the following section we base such qualitative statements deduced from the 
inspection of the distributions (shown as examples in Figs. 8 and 9) on quantitative 
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results of statistical analyses of multivariate distributions, considering also the more 
complicated case of three mass classes (protons: H, carbon: C, iron: Fe). 



4. Non-parametric analyses of the sensitivity to the primary mass 

Non-parametric statistical methods are most efficient and unbiased tools for the analysis 
of multidimensional observable-distributions in order to associate single observed events 
to different classes (say to different masses of the EAS primaries) by comparing 
the considered events with the model distributions without using any pre-chosen 
parameterisation. The methods of decision making and the application to cosmic ray 
data analyses are generally described in Refs. [0, [56|]. They have been outlined for 



studies of muon arrival time distributions in Refs. |27|,|29|. The procedures take into 
account the effects of the natural EAS fluctuations in a quite natural way and are able to 
specify the uncertainties, by an estimate of the true-classification and misclassification 
probabilities. The classification probabilities are determined by the extent to which the 
likelihood functions of the single classes, derived from the simulations, are overlapping. 
Basically the results of such pattern recognition methods, using trained neural networks 
or Bayes decision rules, are dependent from the particular hadronic interaction model 
generating the reference pattern for the data to be studied. In the present analysis 
of the sensitivity of various correlation distributions, prepared from Monte Carlo 
simulations, the so-called one- leave-out test (see Ref. |4"2fl ) is applied, which determines 
the probability that a multidimensional event, taken from the considered (simulated) 
distribution, will be correctly ("true") or incorrectly ("false") classified by the procedure. 



The following figures present the classification (H — > H; C — > C; Fe — ► Fe) 
and misclassification (H — ► C, Fe; C — > H, Fe; Fe —+ H, C) probabilities for various 
correlations of different EAS observable. The discrimination power may be quantified 
by the "Bayes risk" defined as 

R — TT £ ( x i) 

Mi ^ v JJ 

where 




for true classification 

1 otherwise 



of the event Xj of the sample Mj. When considering the results based on simulations, 
one should keep in mind that the primary energy Eq is a-priori known, while in the 
real cases Eq has to be simultaneously estimated from the data, in particular from the 
N e — Nj[ correlation. Within a limited energy range Nj[ can be used as approximate 
energy identifier. From this reason correlations of Nf[ with time observable are not 
specially scrutinised. 

Even with the primary energy given, the fluctuations of N e and N c h or of the 
correlations Ar .50* (-R^) — N e and Ar Q 9 ^Q b (R tJ ) — N ch are so large, that the true 
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Figure 10. Comparison of the classification probabilities deduced from the 
Ar ff 5 ^(i^)-iV e , the Ar * 5 s (ofe (^)-iV e and the Ar^iR^-p^R^-N. correlations 
at the incident energy of (1.00-1. 78)-10 16 eV at = 175-185 m (left). Similarly the 
classification probabilities deduced from the At 9 ^' b (R»)-N ch , the ArU^iR^-Nch 
and the Ar 3 ^ b — p*(i?^) — N ch correlations at the incident energy of (1.00- 
1.78)-10 17 eV are shown (right). The lines are drawn for guiding the eyes. 



classification probability remains modest. Fig. 10 confirms the approximate equivalence 
of the parameter Atq 50 (R^) with the Aro.5o(-R At ) — P* t (^) correlation, though the 
A r o.5o(-R/ t ) — Pu(Rfi) correlation appears obviously preferable in view of the mass 
classification. It is obvious that the quantity p*(Rp) (simultaneously observed with 
multiplicity in arrival time measurements of KASCADE (see Ref. |]34|])) leads to a 
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Figure 11. The influence of the global time parameter AT a 9 ^ b (R^) and the 
accompanying p* on the true and false - classifications of the N e — N^ correlation 
shown for the primary energy (1.00-1.78)T0 16 eV at — 175-185 m. Similarly the 
influence on the true and false - classifications of the N ch — p^R^) correlation is shown 
for the primary energy (1.00-1.78)-10 17 eV (right). For the explanation of the symbols 
see Fig. 10. 



considerable stabilisation of the discrimination effect. Especially when observables (like 
iV e and Ar q ) of comparatively weak discrimination power are correlated. We conclude 
that the time quantities Ar q should be preferentially introduced in combination with 
if ever experimentally possible. 
The effect of global muon arrival time distributions, represented by Ar Q 9 ^'(R tl ) 
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Figure 12. (a) The influence of the global time parameter Atq ^ ob (R^) on the true 
and false - classifications of the N e — iV* r correlation shown for two distances for 
the primary energy range (1.0-1.78)-10 16 eV. For the explanation of the symbols see 
Fig. 10. 



is quantified by the classification probabilities given in Fig. 11. When the detector 
installation is able to measure the N e - Nj[ correlation, adding p*(-R M ) and Ar $f (i^) 
has only a minor discrimination effect, even at larger distances from the shower core. 
This is in contrast to the case, when only N c h and a partial muon number (in our analysis 
approximated by p^R^)) are experimentally accessible. At larger distances, especially 
for the higher primary energies, the correlation of muon arrival times has a clear effect 
of improving the mass discrimination. It is interesting to note that this result holds also 
approximately for local time parameters. 

Fig. 12a displays the independence of the influence of the global time parameter 
Atq J / Q° b (Rn) on the true and false - classifications of the N e — Nff correlation for the 
primary energy range (1.0-1.78)-10 16 eV. At a first glance surprisingly, there is no 
significant tendency with increasing R^. This feature can be understood with the above 
arguments of the minor contribution of the muon arrival time parameter when the (R^ 
independent) N e — N f J correlation can be determined. 

However, due to the generally weaker discrimination power of the N c h — p^R^) 
correlation the influence of Ar 9 ^ b (R IJi ) appears significant in this case (Fig. 12b), 
though also not significantly varying with the distance. For example the improvement 
in terms of the Bayes risk for the case of E = 10 16 eV is only 2%-3% while in the case 
of E = 10 17 eV it is 14%-15%, i.e. not negligible. 

As indicated in Fig. 13 differences in the discrimination effects of different quartiles 
Ar q are hardly obvious and may be obscured by the inherent fluctuations. Studies of the 
classification and misclassification probabilities of EAS observable correlations including 
the age parameters show a non-negligible sensitivity of these EAS parameters. However, 
when alone correlated with N e (N ch ) such a correlation (say N e — s) proves to be a very 
uncertain and instable discriminator due to the fact that fluctuations of both related 
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(1.00-1.78) -10 eV 
175 m< R„ < 185 m N c h-Pn 295 m < R < 305 m 




H C Fe H C Fe 



Figure 12. (6) The influence of the global time parameter Atq J ° (i?^) on the true 
and false - classifications of the N ch — p /1 (R ll ) correlation shown for two distances i? M 
and for the primary energy range (1.0-1.78)-10 17 eV. For the explanation of the symbols 
see Fig. 10. 

observables do obviously corroborate. Adding an observable of another EAS component 
stabilises the classification result. 

From the classification and misclassification studies following features and 
tendencies can be recognised: 

• The reduced time parameter At* = Ar g (i? /1 )/p*(-R^) (where p* represents the 
density of muons with energies > 2.4 GeV) absorbs partially the Ar q (R^) — 
p*(-R M ) correlation and enhances the mass classification sensitivity. 

• The most important features are condensed in Figs. 11 and 12, revealing a 
significant contribution of the muon arrival time information when correlated with 
the EAS observable N c h and p M (i?^), experimentally accessible by the foreseen 
layout of KASCADE-Grande. 

• The shower age parameters s and s ch (which have been scrutinised as additional 
parameters) have modest, but non- negligible effects for the mass discrimination. 

• Correlations of largely fluctuating quantities, resulting from related observations 
(like N e and s, e.g.), lead to a bad classification with large uncertainties, but adding 
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(1.00-1.78) 10 eV 175 m<R< 185 m 




Figure 13. Comparison of the influence of different quartilcs of global arrival time 
distributions on the N c h — correlation. For the explanation of the symbols see 
Fig. 10. 



the correlation with an observable from another EAS component (like At*), even 
also fluctuating, induces an considerable improvement. 

As already mentioned - as additional calculations show - the model dependence of 
the results of the classification probabilities is rather marginal. However, it should 
be remarked, this feature could be in detail different for different EAS parameter 
correlations and is also not yet explored for higher energy ranges. 

5. The efficiency of the observation conditions and the energy and mass 
dependence 

By the observations of muon arrival time distributions special samples of all observed 
EAS events are selected due to the particular observation conditions, especially by the 
energy threshold E thr (= 2.4 GeV in the actual case) of the detected muons and the 
multiplicity threshold n t hr, i-e. the number of muons necessary to define an event with an 
observable arrival time distribution. Thus the event selections are affected by the lateral 
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Figure 14. The probabilities to register EAS events from proton and Fe induced EAS 
of different primary energies under the observation conditions [Ethr =2.4 GeV, n t hr, 
Rfj,) with a 10 • 10 m 2 ideal detector. The variation with i? M and with the multiplicity 
threshold rithr is displayed. 



distribution and the energy spectrum of EAS muon. Consequently the probability being 
accepted in the event sample depends on the observation distance R^, on n thr and on the 
mass of the EAS primary. As discussed in Ref. [34] this feature implies in experimental 
observations of EAS samples efficiency corrections |K| and enables consistency tests by 
varying and nthr i n the measurements. 

Fig. 14 displays the results of calculations of the efficiency to register an EAS event 
under the specified conditions of muon arrival time measurements. 

The efficiency depends strongly on R^ and nthr at lower primary energies (10 15 eV) 
in a qualitatively understandable way. The dependence gets weaker at higher primary 
energies. In general the dependence and the variation with R^ and n^r are different for 
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different primary masses. 
6. Conclusions 

In the present paper based on realistic EAS Monte Carlo simulations the role of 
muon arrival time distributions has been scrutinised in view of correlations with the 
main EAS observables to be measured with the KAS CADE- Grande layout for the 
discrimination of the primary mass at higher primary energies. The main observables 
for that purpose are the total number of the charged particles N c h and a part of 
the EAS muon intensity N^ l art {R fJL ) registered with the KASCADE array, embedded 
in KASCADE-Grande. The partial muon number N^ art {R^) is dependent from the 
location of the EAS centre. For the present simulation studies it has been replaced by 
the density {E^ > 240 MeV). In summary, the analyses of the classification and 

misclassification probabilities, determined with non-parametric statistical techniques 
give evidence that the information from muon arrival time distributions has a potential 
for improving the mass discrimination. This conclusion, being in contrast to findings at 
lower energies, holds for global time quantities as well as for the local ones. Actually 
the EAS profiles of the two different types of time parameters differ by an offset being 
only marginally dependent from the mass. This is due to an increasing "flatness" of 
the < Ati > profile with increasing primary energy, and differences are moved to larger 
distances out of experimental reach. This feature is important with respect to the 
experimental difficulties to define r c and to measure global time quantities. Obviously 
in the considered i? M range the local times provide most of the accessible information. 

The measuring techniques of muon arrival time distributions with the Central 
Detector of KASCADE imply the simultaneous determination of the multiplicity n 
of muons spanning the arrival distributions. This multiplicity is related to the density 
Puffin) °f the higher energy muons detected with the facilities of the Central Detector. 
The quantity p*(R p ) (or the number of muons experimentally detected with the timing 
facility), when correlated with timing measurements (or used as At* = Ar g (i? /J /p* (Rp)) 
plays an important role in improving the effects of the arrival time observations. 

Additionally (as outlined in Ref. p4p as consequence of the specific observation 
conditions muon arrival time observations establish a selected subset of all showers, 



and the determination of the mass composition needs a corresponding correction |60 
Thus the results inferred with a variation of the multiplicity threshold n t h r and of the 
observation distance _R M provide a consistency test for the model used for the non- 
parametric analyses of the observed EAS. 

It should be noted that the present results refer to ideal cases since any detector 
response functions, detector efficiencies and limitations of the statistical accuracies due 
to limited detector sizes have been ignored. Additionally the necessarily limited number 
of simulated EAS implies a limit on the reasonable number of observables used in a single 
analyses. 
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